Localized breathing oscillations of Bose-Einstein condensates in periodic traps. 
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We demonstrate the existence of localized oscillatory breathers for quasi-one-dimensional Bose- 
Einstein condensates confined in periodic potentials. The breathing behavior corresponds to 
position-oscillations of individual condensates about the minima of the potential lattice. We deduce 
the structural stability of the localized oscillations from the construction. The stability is confirmed 
numerically for perturbations to the initial state of the condensate, to the potential trap, as well 
as for external noise. We also construct periodic and chaotic extended oscillations for the chain of 
condensates. All our findings are verified by direct numerical integration of the Gross-Pitaevskii 
equation in one dimension. 

PACS numbers: 05.45.-a, 03.75.Fi, 52.35.Mw, 63.20.Pw 



New techniques for generating Bose-Einstein conden- 
sates (BECs) have opened the door to the investigation of 
a wide range of phenomena and development of concrete 
applications such as atomic interferometers and atom 
lasers. Recent focus has been on BECs trapped in pe- 
riodic magnetic/optical traps the so-called opti- 
cal lattice The addition of a periodic potential to 
the BEC opens the possibility to study exciting dynam- 
ical phenomena of BECs grown in optical traps, such 
as Bloch oscillations and the Josephson effect [|l|,|[D. 
One approach for the study of chains of weakly cou- 
pled condensates, in the tight binding approximation, 
reduces the dynamics of the BECs to the discrete Non- 
linear Schrodinger equation (DNLS) Within 
this approach, each BEC wave function, centered at the 
minima of the potential, is modeled by the same sta- 
tionary wave function with a time dependent amplitude. 
This approximation leads to the DNLS for the local time 
dependent amplitudes. The dynamics described by this 
approach allows for variations in the local population of 
atoms, but are limited to a fixed local distribution of 
the condensed atoms. In this manuscript we use a more 
general approach where the local wave function evolves 
within an ansatz family, allowing for a wider range of dy- 
namical behavior. Using the treatment presented here it 
is possible to include variations in the local populations 
as well as in the shape of the local density: height, width, 
position, phase, velocity, etc. In particular, we focus our 
attention on position-oscillations of each local BEC den- 
sity about its potential minimum by considering a con- 
stant population at each trough of the periodic potential. 
We model the coupled chain of oscillating condensates by 
a nonlinear lattice (a Toda lattice with on-site effective 
potentials) for the oscillation amplitudes. We establish 
the existence of global and localized oscillations of the 
coupled condensates, and demonstrate the ability of con- 



densates to store macroscopic information in localized 
regions of the periodic trap. 

The wave function of a dilute BECs at low temperature 
is governed by the Gross-Pitaevskii equation (CPE) M 
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where 4' = the condensate wave function nor- 
malized with respect to the total number of condensed 
atoms J\f. The atom-atom interactions (considered only 
binary due to the low temperature assumption) are ac- 
counted for by the coupling constant ||] 
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where a is the s-wave scattering length and m the atomic 
mass. The external potential Vcxti'r) is given by the sum 
of the confining potential Vconf (j^) and a periodic optical 
potential Vpcr{r). 

Vc^tir) = Vconf(?^) + Vpci-{r). 




FIG. 1. External potential 14xt with confining (Vconf) and 
periodic (Vpcr(a;)) properties. 

We consider confining traps that give rise to the so-called 
cigar-shaped condensates whose longitudinal direction is 
much larger than their transverse dimension, which is 
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of the order of the heahng length. In this cigar-shaped 
limit it is possible to reduce the GPE equation (|^) to its 
one-dimensional analogue [ic| ] . Our study will be limited 
to the one-dimensional model of attractive BECs. Exten- 
sion of our results to the repulsive BEG in one and higher 
dimensions will be reported elsewhere. Our model applies 
to the dynamics of a chain of quasi-identical condensates, 
in an infinite confining trap. In practice, the harmonic 
external trap has a flat central portion supporting sev- 
eral hundreds of periods of the periodic potential, which 
is sufficient to satisfy the desired near-periodic potential 
(see fig |l|) . It is also possible to load the condensate onto 
T4xt and then adiabatically remove Vconi, leaving only 
the periodic potential. 

After rescaling, the BEG trapped inside a periodic po- 
tential V is governed by the nonlinear Schrodinger equa- 
tion (NLSE): 
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Here corresponds to the rescaled density of 

the condensate, V{x) is the rescaled potential and the 
± sign corresponds to attractive and repulsive (— ) 
atoms. The nondimensionalization is obtained by scaling 
t t/n, X {h/mny/^ x and {mn/4:TTh\a\) \u\^ 

where m is the atomic mass, is the frequency of the 
periodic trap and a is the s-wave scattering length 
We consider a periodic potential V{x), with well spaced 
troughs, constructed by 'concatenating' single well po- 
tentials Vj{x) — Vq tanh^(x — Coy"): 
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V{x) - V,ix)-2NVo. 



The potentials are located at ^q-j — j R where R is 
the potential well spacing. The specific shape of the 
potentials Vj is not central to the construction. How- 
ever, our approximations require that neighboring con- 
densates, grown in adjacent wells, have small overlap. 
This regime corresponds to large enough i? or to a small 
number of atoms per potential well. The large R limit 
can be achieved by tuning the wavenumber of the optical 
trap. Alternatively, new advances in the manipulation 
of condensates permit the construction of traps with ar- 
bitrary spatial patterns by imprinting the pattern in a 
crystal by a grating method or by lithography | |T^ . 

The present manuscript addresses the dynamics of a 
chain of coupled attractive condensates, ie. bright soli- 
tons. For repulsive atoms — dark solitons |ll3|] — a simi- 
lar approach is possible and the detailed analysis will be 
carried elsewhere. Current typical lifetimes for attractive 
condensates range from hundreds to tenths of a second 
before collapse. These lifetimes correspond, for typical 
experimental parameters in *^Rb [Q, to 10'^ rescaled 
time units, long enough to observe 10-20 oscillation pe- 
riods. Moreover, new microtrap techniques and hol- 
low lasers traps |l3] hold out the possibility to increase 



considerably the lifetime of attractive condensates p6[ , 
permitting time to prepare the train of condensates with 
the desired phase configuration. It is also important to 
note that collapse for attractive condensates can be ar- 
rested by controlling the number of atoms. Indeed, if 
the number of atoms in each potential trough is smaller 
than the critical number Nci- associated with the local 
potential, collapse is avoided |]l7|,p|. 

The steady state for a single condensate inside a single 
well potential Vj is a NLSE soliton solution of the form 



We consider the dynamics of a train of condensates 
u{x,t) = '^UQ-^j(x,t) centered at the potential troughs 
^0;j- For simplicity of presentation, we approximate the 
overall dynamics, focusing on two major contributions: 
a) internal dynamics due to interactions between itoy and 
Vj and b) tail-tail interactions between consecutive soli- 
tons. We address each of these terms individually and 
use the linear superposition afforded by soliton pertur- 
bation theory p9| . Within our approximations (small 
tail-tail overlap) it is sensible to discard interactions be- 
tween uoy and Vk for j ^ k and interactions between 
non-nearest neighbor condensates |19|. 

The evolution, for small perturbations, of the steady 
state soliton uq-j inside its on-site potential Vj can be 
approximated by an oscillating soliton ansatz uj(x, t) = 
UQ-j {x — , t) with constant height and width |20(| . The 
position ^j (t) for each condensate is then described by a 
particle inside an effective potential: 



= -Kff (0 - Co;,)- 



(4) 



The effective potential VcS may be obtained by soliton 
perturbation techniques or alternatively by demanding 
that the evolution of ^j (t) in (||) respects the invariance 
of the the NLSE (|) energy 
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It is well known pO] that the effective potential 
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is proportional to the overlapping integral between the 
displaced BEG density (^) and the on-site potential. 
Note that in (^ the on-site potential Vj is centered at 
X — ^o;j and the BEG density \uQ-j{x — C)\^ is centered 
a.t X — ^o-j + ^ (see (||)). The exact form of Vcs de- 
pends upon the on-site potential. For the particular 
case under consideration, Vj{x) = Vq tanh^(a; — ^o-.j), 
the effective potential admits the expansion Vcs{x) oc 
j^x'^ — -^x^ + o{x^). More generally, the effective poten- 
tial can be approximated, for small oscillations, by 
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where a and [3 encode the shape information of Vj. In 
particular, even symmetry of Vcs is inherited from Vj. 

Tail-tail interactions for neighboring condensates re- 
sult in a complex version of the Toda lattice involving 
position and phases ||2l|. It is possible to further reduce 
the dynamics under the assumption that relative phases 
for consecutive condensates are constant. This is a rea- 
sonable approximation in the regime where the solitons 
are kept well apart p^ ]. In particular we consider the 
case of a TT phase shift between consecutive condensates, 
since in the steady-state (^j = 0) this configuration re- 
duces to the Jacobi elliptic cosine which is known to be 
stable |18|. In contrast, a chain of condensates with zero 
phase shift reduces to the third Jacobi elliptic function 
(dn), which is unstable |Q. With a relative phase of tt 
the tail-tail interactions reduce to a real Toda lattice (TL) 
Q on the positions: '(^ = 4 (e~(^^~«^-i) - e~(^^+i~^^)). 
In practice, the tt phase shift can be implemented by 
phase design on the initial configuration of the conden- 
sates . From now on we consider that the amplitudes 
of oscillation of the condensates are small. Thus, since 
the steady state is stable, we eliminate the possibility of 
a condensate hopping to a neighboring lattice trough. 
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FIG. 2. Left: Homoclinic tangle for the two-dimensional 
map (0). The homoclinic orbit {..., P-i, Po, Pi, ...} belongs 
to the stable (sohd line) and unstable (dash line) manifolds. 
Right: The corresponding configuration for the oscillation 
amplitudes. {R = 10, w = 17.671 and Vo = 0.1.) 

Exploiting the linearity of soliton perturbation theory 
we combine the on-site potential <Q) with the tail-tail 
interactions to find a lattice differential equation (LDE) 
on the condensate positions 

- 4 (e-fe-?.-) - - v:^{^, - eo,,). (8) 

To find oscillatory solutions to ^ we use an oscillat- 
ing ansatz p5[ |, writing the position of the jth con- 
densate as a combination of oscillatory modes £,j{t) = 
J2k^ji^) cos{kujt) centered at Aj{0) — ^oj- For small 
oscillations a one-mode ansatz 



(t) = 6o.j + Aj cos(wi) 



(9) 



is sufficient to capture the essential dynamics. We sub- 
stitute (^ into (^), expand and match terms to find a re- 
currence relation between the oscillation amplitudes Aj: 



hA^A, 
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(10) 



with a = 2— ti)^-f a e^' j ^ and b — 3 (3 /IQ. We introduce 
Hn = An and Xn = An-i and recast the second order 
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recurrence relation {nw as the 2D map M: 



M 



yn+1 = {a + byl)yn- x„ 



Up to the approximations made to obtain equation (|1C 
solutions of (|l^) prescribe the amplitudes representing 
oscillatory solutions for the condensates. We remark that 
the oscillation frequency uj and the form of the on-site 
potential are incorporated through the parameters a and 
6of (|l^). In particular, one expects families of oscillatory 
solutions parameterized by their frequency uj. 




FIG. 3. Localized breathing oscillation in a chain of weakly 
coupled condensates in a periodic potential. This local- 
ized oscillation is obtained by full numerical solution of the 
Gross-Pitaevskii equation with an initial condition prescribed 
by our dynamical reduction. Only condensates with index 
through 4 are shown (the oscillations are symmetric with 
respect to the condensate with index 0). The bottom plane 
depicts ut — darker areas corresponds to regions where u{x, t) 
undergoes greater temporal variation. (_R = 9, a; « 0.11388 
and Vo = 0.025.) 

We are interested in constructing localized breathers 
(localized oscillations) for the condensate dynamics. 
These solutions correspond to orbits homoclinic to the 
origin for the recurrence map (pT|). The correspond- 
ing orbits satisfy Aq ^ and lim„_>-|-oo An — with 
an exponential decay rate. A necessary condition for 
the existence of a homoclinic orbit is hyperbolicity of 
the origin. This is satisfied for all |a| > 2 or equiva- 
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> 2. This condition provides 
constraints on the inter-site spacing R and breather fre- 
quency UJ for a given on-site trap Vj, which determines 
a. We find numerically, for reasonable values of i? (i? « 
a few condensate widths), that the stable (PF*) and un- 
stable manifolds (VF") of the origin intersect in a ho- 
moclinic tangle (see left panel in Fig. ||). Consider the 
trajectory of the two-dimensional map M that starts at 
the intersection point Pq G n IF" (see left panel in 
Fig. H). From the invariance of the stable and unsta- 
ble manifolds, each forward and backward iteration of 
Po, labeled {. . . , P_2, P-i, Po, Pi, P2, ■ • ■} on Fig. |, lies 
in the intersection. From area orientation it follows that 
this trajectory takes each second intersection pq]. The 
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amplitudes for the localized oscillations are then given 
by the ordinates (j/„ = An) of the homoclinic orbit (see 
right panel in Fig. ||). 

The homoclinic orbit of the 2D map induces a breather 
on the condensate positions through the ansatz (j^). An 
example of such localized oscillations for the condensates 
is depicted in Fig. ^, in which a central condensate oscil- 
lates with a maximal amplitude Aq and the condensates 
on either side oscillate with amplitudes A±n which de- 
crease exponentially in increasing n (see right panel in 
Fig. 0). The asymptotic decay rate for the oscillations 



(AL^"^' = A^'"') is prescribed by the eigenvalues at the ori- 
gin A± = (a ± Va^ — 4)/2. Note that, because ( |Tl| ) has 
a; <-*■?/ symmetry, the breather is symmetric with respect 
to the central condensate (A_ = •^+^)- We stress that 
the solution depicted in figure ^ is obtained by numerical 
integration of the full Gross-Pitaevskii equation (|l|) from 
initial conditions prescribed by the homoclinic orbit of 
our reduced 2D map. Due to the approximations used in 
our approach, orbits of the reduced 2D map are not ex- 
act solutions of the full GPE equation(|^). Nevertheless, 
if we restrict our attention to small oscillations for weakly 
coupled BECs, there is a good correspondence between 
the reduced dynamics and the original Gross-Pitaevskii 
equation. This correspondence is reinforced by the struc- 
tural stability of the orbits for the reduced 2D map (see 
below) . 
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FIG. 4. Robustness of the perturbed localized breather in 
the Gross-Pitaevskii equation: a) unperturbed and b) per- 
turbed, ai and bi depict a density plot of \ut{x,t)\ — darker 
areas correspond to larger u{x, i) variations. Insets a.2 and 
as depict Ut{x,to) and ln(|itt(a;, to)|) ai t = to for the un- 
perturbed system, whereas insets b2 and ba depict the cor- 
responding quantities at ti for the perturbed systems. Same 
parameter values as in Fig. |^. 

Localized breathers represent an important structure 
in the dynamics of the GPE equation describing a mech- 
anism whereby energy and information can be pinned 
down on the periodic potential lattice. The occurrence 
of localized breathers in weakly coupled oscillatory units 
is a common phenomenon. Indeed, the existence of lo- 
calized breathers has been formally established in gen- 
eral nonlinear Hamiltonian lattices of weakly interacting 



oscillators [|2^ by continuation from the so called an- 
ticontinuum limit corresponding to the uncoupled case 
(i? = +oo) The structural stability of the homo- 

clinic tangle, arising from the transverse nature of the 
intersections of the stable and unstable manifolds, im- 
plies the persistence of the breather solution of the 2D 
map ( pd] ) under parameter changes, in particular insur- 
ing the existence of the breather solution for the GPE. 
This persistence, together with the dynamic stability of 
the steady state solution {An = 0), ensures the existence 
of the breather solution for the GPE. Therefore, despite 
the various approximations used in the reduction of the 
GPE to the 2D map (pd]), we observe localized breather 
oscillations in the original GPE dynamics (Fig. ||) for 
the predicted parameter values. More surprising is the 
robustness of the breathing dynamics to significant per- 
turbations under the GPE dynamics (Fig. ^j. 

To demonstrate this robustness we construct the ini- 
tial configuration (IC) as a concatenation of solitons with 
velocities, heights, widths and positions predicted by our 
analysis, and with a tt phase shift between consecutive 
solitons (as in Fig. ^). We perturb each soliton in the 
chain by adding a random value to each of the initial 
velocities, heights, widths, positions and phases. The 
bounds for the perturbations correspond to 

a) 15% of the maximal velocity on velocities, 

b) 10% on heights, 

c) 10% on widths, 

d) 15% of the maximal amplitude (0.15 xAq) on posi- 
tions and 

e) 15% of TT on phases. 

Note that the perturbations on positions and velocities 
are proportional to the collective maximum and not to 
the individual values. After adding the perturbations we 
concatenate the solitons and integrate the NLS (||) using 
a pseudo-spectral method. Additionally, we include 

f) 5% (stationary) perturbation to the potential and 

g) a noise level of 10^^ to u{x, t) at every time step of 
the integration. 

The total run is about 2 x 10^ iterations. Fig. depicts 
the evolution of Ut{x,t) for the unperturbed (left) and 
the perturbed (right) breather from direct simulations 
of the GPE equation. Despite the large perturbation to 
the IC and the additive noise, after a brief transient, the 
breather settles down and retains its localization with 
an approximate exponential decay (inset bs). It is inter- 
esting to note that stronger perturbations to u did not 
necessarily destroy the localization but often resulted in 
a slow excursion of the localized region along the lattice. 
The possibility of breather mobility in nonlinear lattices 
triggered by perturbations has been investigated previ- 
ously [pi. 
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The robustness of the locahzed breather dynamics pre- 
sented above opens the possibihty for experimental cor- 
roboration. The large perturbations to the IC result in an 
initial transient involving radiative losses, which would 
correspond experimentally to a small number of atoms 
being spilled away from the central cloud and absorbed 
by the non-condensed atoms at the periphery. After the 
transient, the breather settles down and retains its local- 
ization with an approximate exponential decay (inset bs). 
While the breathers we constructed are robust to a wide 
host of perturbations, breathers are non-generic in the 
sense that arbitrary initial conditions do not necessarily 
produce a breather. Also, with large enough perturba- 
tions the breathing phenomena are entirely destroyed. 
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FIG. 5. Typical phase space for the reduced 2D map (pT[). 
Left: some periodic orbits and the homoclinic tangle. F and 
F' are fixed points, {Gi, G2} is a period 2 orbit. Some higher 
order periodic orbits: period 3 (triangles and filled squares) 
and period 4 (open squares). Right: behavior near the fixed 
point F (diamond). The map displays quasiperiodic orbits 
which disappear at the separatrix originating from the mani- 
folds of a period 11 orbit (crosses). Outside this separatrix a 
single chaotic orbit is depicted (small dots). 

We may use the dynamical reduction described above 
to devise other breathing phenomena of Within our 
approximations, any bounded nontrivial orbit of the 2D 
map M (|l^) gives rise to complex oscillatory behav- 
ior of the condensates in (^). In particular, periodic 
points of M correspond to global oscillations. The map 
M has three fixed points: the origin, F = {x*,x*) and 
F' = {-x*,-x*) (see Fig. |), where x* = y^{2- a)/b. 
The fixed point at the origin gives rise to the trivial 
stationary solution A^s— 0. The fixed points F and 
F' yield solutions in which all condensates oscillate in 
phase with the same amplitude x* . The correspond- 
ing global breather for the original system (||) is de- 
picted in Fig. ^.a. Other interesting orbits arise from 
higher order periodic points. For example, the period 
two orbit G2 = -M(Gi) = ]VP{G2) (see F ig. |) where 
Gi = {x,—x) and G2 = {—x,x) with x = -y^— (2 -I- a)/b. 
This period two orbit corresponds to an amplitude config- 
uration {. . . , X, —X, X, —X, . . .}, ie. anti-phase oscillations 
(see Fig. y.b). It is in principle possible to construct more 
complex patterns for the global oscillations from higher 
order periodic orbits of the 2D map (cf. Fig. ||). 

The 2D map also predicts global oscillations which are 
quasi-periodic in site index n. These orbits exist near 



the fixed point F (see Fig. ||, right panel). The corre- 
sponding breather for the full periodic NLSE (||) has an 
in-phase global oscillation with a small modulation of the 
amplitudes (given by the rotation number of the quasi 
periodic orbit around F). An interesting possibility for 
the dynamics of coupled condensates is the prospect of 
chaotic evolution. In a neighborhood of the fixed point 
F there is a region containing chaotic orbits correspond- 
ing to chaotic oscillations for the condensates (see right 
panel in Fig. ||). This corresponds to chaotic oscillations 
for the condensates (see Fig. |^.c). It should be possible, 
in principle, to find the onset of chaos for the conden- 
sates by analyzing in more detail the reduced dynamics 
1h 
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FIG. 6. Global oscillations of the condensates from sim- 
ulations of the GPE equation. These global oscillations 
arise from periodic points of the reduced dynamics (^l|). a) 
In-phase oscillations corresponding to the nontrivial fixed 
point F (see Fig. ^). b) Anti-phase oscillations corresponding 
to a period two cycle of (llll) with alternating amplitude sign. 
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FIG. 7. Chaotic global oscillations of the condensates from 
GPE dynamics corresponding to the chaotic orbit of the re- 
duced map (|ll|) (see Fig. P). 

We have constructed a variety of novel global and local- 
ized oscillatory behaviors of BECs in periodic potentials, 
identifying these solutions with orbits of a reduced 2D 
map. A key ingredient of the construction of localized 
oscillations is the existence of a homoclinic tangle. We 
demonstrate the surprising robustness of these solutions 
to perturbations. Since BEG experiments are quite del- 
icate, we do not expect that direct manipulation could 
produce an exact initial condition corresponding to a lo- 
calized oscillation. Nonetheless, we believe that localized 
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oscillations may be observed in weakly coupled conden- 
sates which are appropriately engineered and then per- 
mitted to radiate away spurious energy. The techniques 
presented here can, in principle, be extended to lattices in 
higher dimensions such as vortex lattices. This has possi- 
ble implications to modeling the interactions of atoms in 
optical traps that could potentially be used for quantum 
computing po|| . 
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